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0\ ■ Abstract 

p || We update the Standard Model predictions for the rare decays tt+z/z/ 

Q_i' and Kl tx^uu. In view of improved hmits on Bg-Bg mixing we derive a strin- 

1^ . gent and theoretically clean Standard Model upper limit on B[K^ —>■ tc'^uu), 

which is based on the ratio of B^-Ba to Bg-Bg mixing, AM^/ AMg, alone. This 
method avoids the large hadronic uncertainties present in the usual analysis of the 
X ; CKM matrix. We find B{K+ — > n^vu) < 1.67 ■ 10 which can be further im- 

■ proved in the future. In addition we consider the extraction of \Vtd\ from a future 

measurement of B{K~^ — > vr+z/z/), discussing the various sources of uncertainties 
involved. We also investigate theoretically clean constraints on B{Kl ir^uu). 
We take the opportunity to review the next-to-leading order (NLO) QCD cor- 
rections to ^ TTuu, Kl fJ'^^f^'', B — s> Xuu and B including a 
small additional term that had been missed in the original publications. The 
phenomeno logical impact of this change is negligible, the corresponding numeri- 
cal shift being essentially within the small perturbative uncertainties at the NLO 
level. 



1 Introduction 



The rare decays ty^uu, Kl it uu, B Xs^d^^v and Bs^d are 

very promising probes of flavour physics. They are sensitive to the quantum 
structure of Standard Model flavour dynamics and can at the same time be com- 
puted theoretically to an exceptionally high degree of precision. In particular the 
kaon modes — > tc^uu and Kl — > tt^uu have received considerable interest in 
recent years. The theoretical status of these decays has been reinforced by the 
calculation of next-to-leading order QCD corrections The corrections to 

the top-quark dominated modes — >■ n^uu, B — >• Xg^d^u and Bg d 0) 

have recently been recomputed in ^j. The results of |0 confirm the calculations 
of ||1], 3 up to a small term arising from the box diagram that had been overlooked 
in [0 . We are in agreement with that this term needs to be included to obtain 
the complete NLO correction. As will be explained in more detail below, the 
missing piece follows from a simple one-loop calculation, which can be performed 
separately, and it can be added to the old result. The latter was obtained from 
the two-loop matching calculation in 0, which remains unaffected. 

In contrast to — > ir^ui?, where only the top-quark contribution is relevant, 
also the charm sector is important for tt^uu. This sector, which has not 

been considered in 0], will also be discussed in the present article. 

The numerical impact of the modification is small. The changes are within 
the perturbative error of the NLO calculation, both for the top and the charm 
contribution. All qualitative conclusions, regarding the NLO error analysis from 
residual scale dependence and the theoretical precision that can be achieved in 
the calculation of these decays, remain unchanged. 

The purpose of this paper is twofold. We will first present the complete NLO 
effective hamiltonians for tt^uu (section 2), the short- distance part of 

Kl — > (section 3) and the top-quark dominated modes Kl ti^uu, B 

Xg^di^u and B^^d ^^^^ (section 4), including the additional NLO contribution 
from the box diagram. The origin and the derivation of this contribution are 
explained in the Appendix. Secondly, we update the Standard Model predictions 
for B{K^ —>■ TT^uu) and B{Kl —>■ tt^uu) and discuss some phenomenological 
aspects of these decays that have not been emphasized before. In particular, we 
derive in section 5 a clean correlation between the ratio of Bd-Bd and Bg-Bg 
mixing and the K~^ —>■ tc'^uu branching fraction. In view of recent experimental 
progress, both for K^ n^uu and B-B mixing, this could provide an interesting 
Standard Model test in the near future. We also analyse the determination of 
\Vtd\ from B{K^ ti^vv) and investigate the theoretical uncertainties involved. 
Section 6 explores the impact of AMd/ AMg and sin 2/5, as obtained from the 
CP asymmetry in Bd{Bd) J/'^Ks, on the neutral mode Kl — >■ n^vv. The 
combination of these three theoretically clean observables offers another stringent 
consistency check on the Standard Model. We conclude in section 7. 



1 



2 Effective Hamiltonian for K 



The effective liamiltonian for vr+z/z/ can be written as 

neff = ^ 2nsi^'ew , ^ {^*sVcAl + V:,Vt,X{x,)) i7sd)v^Aim)v-A . (1) 

Tfie index l=e, /i, r denotes tlie lepton flavour. Tlie dependence on tlie cliarged 
lepton mass, resulting from tlie box-graph, is negligible for the top contribution. 
In the charm sector this is the case only for the electron and the muon, but not 
for the r-lepton. 

The function X{x), relevant for the top part, reads to 0{as) and to all orders in 
X = m'^/M^ 

(2) 
(3) 



a. 



X{x)=Xo{x) + -^Xi{x) 
An 
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where x^ = fi^/M^ with = 0{mt) and 
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XI 
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(4) 



(5) 



The //-dependence in the last term in cancels to the order considered the 
//-dependence of the leading term Xo(x(/i)). 

The expression corresponding to X{xt) in the charm sector is the function X^^. 
It results from the RG calculation in NLLA and is given as follows: 



Xnl — Cnl 



AB 



(1/2) 
NL 



(6) 



Cnl and B^j^^ correspond to the Z°-penguin and the box- type contribution, 
respectively. One has 
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Here K2 = K^^^'^^, m = nic, r = mf/ml(fi) and m; is the lepton mass. We will 
at times omit the index / of X^^. In (^ - (0) the scale is = C(mc). For the 
charm contribution we need the two-loop expression for as{fi) given by 
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(12) 



The effective number of flavours is / = 4 in the expressions above for the charm 
sector. The QCD scale in ([n|) is A = = (325 ± 80) MeV. 

Again - to the considered order - the explicit ln(/i^/m^) terms in (^) and (p!0D 
cancel the /i-dependence of the leading terms. 

For phenomenological applications it is useful to define for the top contribution 
a QCD correction factor r]x by 



X{xt)=rix-Xoixt) Vx = 0.99A 



(13) 



With the MS definition of the top-quark mass, rrit = mt{mt), the QCD factor rjx 
is practically independent of mj. This is the definition of rrit that we will employ 
throughout this paper. 

For the charm sector it is useful to introduce 



Po{X) = ^ 



(14) 



with the Wolfenstein parameter P] A = 0.22. Results for the charm functions are 
collected in tables |1] and 0, where rric stands for the MS mass mc{mc)- 



3 Effective Hamiltonian for (Kl )sd 

The analysis of {Kl — > f^^fJ'~)sD proceeds in essentially the same manner as for 
TT^vv. The only difference is introduced through the reversed lepton line 
in the box contribution. In particular there is no lepton mass dependence, since 
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Table 1: The functions and X^-^ for various Aj-^ and m,,. 







4 




-4 


A^[MeV]\m, [GeV] 


1.25 


1.30 


1.35 


1.25 


1.30 


1.35 


245 


10.73 


11.60 


12.50 


7.34 


8.06 


8.81 


285 


10.44 


11.31 


12.20 


7.05 


7.77 


8.52 


325 


10.14 


11.00 


11.90 


6.75 


7.47 


8.21 


365 


9.82 


10.69 


11.58 


6.44 


7.15 


7.89 


405 


9.49 


10.35 


11.24 


6.10 


6.81 


7.55 



Table 2: The function Po{X) for various and rric. 







A^[MeV]\mc [GeV] 


1.25 


1.30 


1.35 


245 


0.410 


0.445 


0.481 


285 


0.397 


0.432 


0.468 


325 


0.385 


0.419 


0.455 


365 


0.371 


0.406 


0.442 


405 


0.357 


0.391 


0.427 



only massless neutrinos appear as virtual leptons in the box diagram. 

The effective hamiltonian in next-to-leading order can be written as follows: 
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V22 

The function Y{x) is given by 
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Similarly to (|T^ one may write 

Y{xt)=VY-Yo{xt) r]y = 1.012 (19) 

The RG expression Y^l representing the charm contribution reads 

Ynl = Cnl - B^j^l'^^ (20) 

where C^vl is the Z°-penguin part given in (|^ and B^^^^"^^ is the box contribution 
in the charm sector, relevant for the case of final state leptons with weak isospin 
T3 = -1/2. One has 



Table 3: The functions Inl and Po{Y) = Inl/A"^ for various and rric. 





>nl/10- 


^4 


PoiY) 


A^[MeV]\m, [GeV] 


1.25 


1.30 


1.35 


1.25 


1.30 


1.35 


245 


2.75 


2.94 


3.13 


0.117 


0.125 


0.134 


285 


2.80 


2.99 


3.19 


0.119 


0.128 


0.136 


325 


2.84 


3.04 


3.24 


0.121 


0.130 


0.138 


365 


2.89 


3.09 


3.29 


0.123 


0.132 


0.140 


405 


2.92 


3.13 


3.33 


0.125 


0.133 


0.142 
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Note the simple relation to B^^^j^^ in ( |lOl) (for r = 0) 



(21) 
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(22) 



4 Effective Hamiltonians for Kl - 

Xs^-uu and Bs4 l^l' 



7r"z/z/, B 



With the above results it is easy to write down also the effective hamiltonians for 
Kl — > tt'^uu, B — s> Xs^^z/z/ and Bg^d l^l^ . Since only the top contribution is 
important in these cases, we have 



n 



eff 



a/2 2tt sin^ Giy 



V*^Vtn'X{xt)inn')v-A{^i^)v-A + h.c. 



(23) 
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for the decays Kl — ^ ti^vv, B Xsvv and B X^i/z/, with {nn') = (sd), (bs), 
(bd), respectively. Similarly 

^eff = -"n V:^Vtn'Yixt)inn')v.Aill)y.A + h.c. (24) 

y z ^7T sm 

for Bg — s> and B^ l^l^, with (nn') = (6s), (bd). The functions X, y are 
given in (^ and (|T^). 

5 Phenomenology of tt+z/z/ 

5.1 General Aspects and Standard Model Prediction 

The branching fraction of K'^ n^uu can be written as follows 



B{K+ TT+Z/Z/) 



2 /t^ ^ , \ 2 



(^^X(..)) +(^Po(X) + ^X(..) 



(25) 



^+ = '2 • 4n = 4-11 • 10-" (26) 

27r"= sm fc)vK 

Here Xt = mf/M^r, Xi = V*.Vid and rx+ = 0.901 summarizes isospin breaking 
corrections in relating —>■ tt^vv to the well measured leading decay — > 
7r°e^z/ 0. In the standard parametrization Ac is real to an accuracy of better 
than 10"'^. We remark that in writing B{K^ tt^uu) in the form of (pS]) a 
negligibly small term ~ (X^^ — ^nl)"^ has been omitted (0.2% effect on the 
branching ratio). 

A prediction for B{K~^ tc'^uu) in the Standard Model can be obtained 
using information on kaon CP violation (ek), \Vub/Vcb\ and B — B mixing to 
constrain the CKM parameters ReA( and ImAj in (p5|). This standard analysis 
of the unitarity triangle is described in more detail in Here we present an 
updated prediction using new input from the 1998 Vancouver conference We 
take 

Ex = 0.80 ±0.15 p^fBsV^ ^ ^ 2 Po(^) = 0.42 ± 0.06 (27) 



for the (scheme-invariant) kaon bag parameter B^, the SU{3) breaking among 
the matrix elements of Bg Bg and B^-B^ mixing, and Po{X) in (125|) , respectively. 
Next we use [|T^ 

mt = Tht{mt) = (166±5)GeV = 0.040±0.003 |Kfe/Kfe| = 0.091±0.016 

(28) 



and the experimental results on Bd — Bd and Bg — Bg mixing |10 



AMd = (0.471 ± 0.016)ps"^ AM, > 12.4ps-^ (29) 
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Scanning all parameters within the above ranges one obtains 



B{K+ TT+z/z/) = (0.82 ± 0.32) ■ 10 



-10 



(30) 



where the error is dominated by the uncertainties in the CKM parameters. Eq. 
(BUI) may be compared with the result from Brookhaven E787 [p!T| 



(4.2 



(31) 



Clearly, within the large uncertainties, this result is compatible with the Standard 
Model expectation. 



5.2 Upper bound on B{K+ tt+uu) from AM^/AMg 

Anticipating improved experimental results, the question arises of how large a 
branching fraction could still be accomodated by the Standard Model. In other 
words, how large would B{K^ 7r+z/z/)exp have to be, in order to unambiguously 
signal New Physics. In this context we recall that the clean nature of —>■ vr+z/z/ 
implies a relation between the branching ratio and CKM parameters with very 
good theoretical accuracy. However, in order to constrain the poorly known 
CKM quantities to predict (^0|) , one introduces theoretical uncertainties (related 
to iKjft/Kbl or Bk) that are not intrinsic to tt'^uu itself. We would 

therefore like to investigate to what extent an upper bound can be derived on 
B[K^ — i> TT^vv), without relying on iKib/Kftl or the constraint from kaon CP 
violation {sk) involving B^- For this purpose we will make use of the ratio 
AM(i/AM<.. This is motivated by the theoretically fairly clean nature of this 
ratio and the improved lower bound on AM^ (|29D. 

In terms of Wolfenstein parameters \, A, q and 7] |H| one has p[ 



Vt 



td 



V, 



ts 



1 + A2(2^- 1) 



(32) 



where 12 




2 ) V 2 ; 

A measurement of AM(i/AMs determines Rt according to [§] 



(33) 



^^-TVamT^ 



2g) 



r.. = W^^^^ (34) 



The ratio of hadronic matrix elements r^^ has been studied in lattice QCD. As 



discussed in [jT3[-||T5[, the current status can be summarized by 



1.14 ±0.08 



(35) 
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This result is based on the quenched approximation, but the related uncertainties 



are expected to be moderate for the ratio (|35D. In the following we shall use 



r,, = 1.2 ± 0.2 (36) 

Since r^^ = 1 in the 5'[/(3)-flavour symmetry limit, it is the difference r^^— 1 that is 
a priori unknown and has to be determined by non-perturbative calculations. To 
be conservative, we have assigned a 100% error on the 5'f/(3)-breaking correction 
in (H). 

We may next rewrite (pSj) using an improved Wolfenstein parametrization |12 
as i 

B{K+ Ti+uu) = K+A^X\xt)- \{cyf]f + (f?o - Q?] (37) 

Eq. (|37D defines an (almost circular) ellipse in the (^, ff) plane, centered at (^o, 0). 
Now, for fixed the maximum possible branching ratio occurs for 77 = and is 
given by 



5(K+ n^iyiy)max = — 

a 



Po{X) + A'X{xt 



A V AM, 



(39) 



from (|32|) , (|3^) and (|37|) . This equation provides a simple and transparent relation 
for the maximal B{K^ — ^ t^^uv) that would still be consistent with a given value 
of AMd/ AMg in the Standard Model. In particular, a lower bound on AMg 
immediately translates into an upper bound for B{K^ — > tt^uu). We stress that 
( p9D is theoretically very clean. All necessary input is known and does not involve 
uncontrolled theoretical uncertainties. At present the largest error comes from 
Tsd, but a systematic improvement is possible within lattice gauge theory. 

For these reasons (^) can serve as a clearcut test of the Standard Model and 
has the potential to indicate the presence of New Physics in a clean manner. 
Using 



V AM 
one finds 



^^'^ < 0.2 A< 0.89 PoiX) < 0.48 X{xt) < 1.57 r.^ < 1.4 



(40) 



B{K+ ^ = 1.67 ■ 10"!° (41) 

This limit could be further strengthened with improved input. (For central values 
of parameters the number 1.67 in ( ^Tj ) would change to 0.94.) However, even with 
present knowledge, represented by (^TJ), the bound is strong enough to indicate 
a clear conflict with the Standard Model if B{K^ tt^uu) should be measured 
at 2 ■ 10~i°. 

In table ^ we illustrate how the upper bound on B{K^ n^uu) (bound A) 
improves when the maximal possible value for Vgd is assumed to be smaller. The 
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outcome is compared with the result from the standard analysis of the unitarity 
triangle (bound B). We find that as long as Vds is higher than 1.4, the r^s bound 
has no impact on the maximal value in the standard analysis. We also observe 
that bound A, which is theoretically very reliable, is only slightly weaker than 
the less clean result of bound B. 

Table 4: B{K+ vr+z/i?)„„^ ■ 10^° from AM^/AM, alone (bound A), compared 
with the maximum value from the standard analysis (bound B), where also the 
information from iKib/Kbl, Ek and AM^ is used. The bounds are shown for 
various rsd,max, the maximum of the SU{3) breaking parameter rsd- 



fsd,max 


Bound A 


Bound B 


1.40 


1.67 


1.32 


1.30 


1.49 


1.21 


1.25 


1.40 


1.17 


1.20 


1.32 


1.14 



5.3 \Vtd\ from B{K+ tt+uu) 

Eventually, a precise experimental determination of B{K^ tt^uu), in particu- 
lar if compatible with Standard Model expectations, can be used to extract \Vtd\ 
directly from (^) 0. We would like to illustrate such an analysis here by 
detailing the sources of uncertainty and their impact on the final result. Our 
findings are summarized in table We remark that the sensitivity of \Vtd\ to 
variations in the input is fairly linear for the parameters B{K^ vr+z/z/), and 
Vcb through A^, so that the effect of other choices for the errors can be easily 
infered from this table. 



6 Phenomenology of Kl 



The rare decay mode Kl — ^ n^uu is a measure of direct CP violation [|T7l and 
therefore of particular interest. Using the effective hamiltonian ( PBD and summing 
over three neutrino flavours one finds |H, B 



Kl = K+- 



rxi^ t{Kl) 



1.80 ■ 10"^° 



(42) 



(43) 



'tk+t{K+) 

with tk^ = 0.944 the isospin breaking correction from and /€+ given in 
Using the improved Wolfenstein parametrization we can rewrite (^2]) as 



B{Kl ^ ttVz/) = KLr]^A'X\xt 



(44) 
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Table 5: Sensitivity of \Vtd\, extracted from vr+z/z/, to the branching ratio, 

to Standard Model parameters, and to the renormalization scales in the top and 
charm sector, fit and fic respectively. The latter measure the residual theoretical 
uncertainties from perturbation theory at NLO. rrit and rric are MS masses. The 
branching ratio is assumed to have the value shown in the table for the purpose of 
illustration. Variation of the inputs within the indicated errors yields the quoted 
percentage change in \Vtd\ around \Vtd\ = 1.029 ■ 10~^, corresponding to central 
parameter values. 





fit/GeV 


/ic/GeV 


(1.0 ±0.1) • 10-1° 


100-300 


1-3 


±6.8% 


±0.5% 


±4.5% 




\Vub/V,b\ 


mt/GeV 


mJGeY 




0.040 ± 0.003 


0.091 ±0.016 


166 ±5 


1.30 ±0.05 


0.325 ±0.080 


±7.6% 


±0.6% 


±3.5% 


±2.9% 


±2.2% 



Using the same procedure (standard analysis) and input as for the derivation of 
(pOD, one obtains 

5(/s:l ^ vrVz/) = (3.1 ± 1.3) • 10"" (45) 

Again the error comes almost entirely from the uncertainties in the CKM param- 
eters. 

As in the case of tt+z/z/ discussed above, one may ask the question 

about the maximal —>■ tc^uu branching ratio that could be tolerated within 
the Standard Model. A value of 4.4 ■ 10"^^ is suggested by (|45|) . Relying on the 
standard analysis of the CKM matrix, this result necessarily involves hadronic 
uncertainties (related for instance to Ek and |Ki6/Kfe|). It would be desirable 
to have a relation that could yield comparable information, but based only on 
quantities with very good theoretical control. A promising possibility is again 
AMd/ AMg, in conjunction with sin 2/3, to be obtained from the CP asymmetry 
in Bd{Bd) J/^Ks- For all practical purposes the latter quantity has no 
theoretical uncertainty, like the Kl tt^uj? branching fraction itself. Using the 
CKM relation 

f = ^ (l - ^1- sin' 2(3^ (46) 

together with (0) gives 

BiK, n'uu) = ^^LA'X\xt)^^^^ (l - /T^^l^) (47) 

Here we have dropped the small (9(A^) term in (|^) for simplicity. It could 
always be included if a still higher accuracy should be required. In any case, the 
omission of this correction is conservative when we are interested in an upper 
limit on B{Kl — > tt^uu) (unless g > 0.5). A useful numerical representation of 
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B{Kl vrVz/) = 12.0 ■ 10"" 



0.04/ V166GeV/ V1.2/ 0.04 



1 - VI -sin^2/5j (48) 

Experimental results for AM^/AMg and sin 2/5 immediately translate into a pre- 
diction for B{Kl — >■ TT^i'i'). If only upper bounds are available for AMd/AMg 
and sin 2/5, still a clean upper bound on B{Kl tc^uu) can be obtained. 

At present, with AM^/AM^ < 0.04 and no information on CP violation in B 
decays (sin 2/? < 1), this bound is not yet very strong as can be seen from (^8]). 
Note however that the last factor in (^7]), (^) decreases rapidly when sin 2/3 is 
restricted to values smaller than unity. Improved results on AMd/AMg will also 
contribute to make this upper bound competitive with the standard analysis (^5]). 
For example, taking Vet < 0.043, nit < 171 GeV, r^a < 1-4, AMd/AM^ < 0.04 
and assuming sin 2/3 < 0.7 yields 

B{Kl Ti^vv) < 6.7 ■ 10-" (49) 



Let us emphasize again that the main virtues of (^Tf ) with respect to the stan- 
dard CKM analysis are the very high reliability of the theoretical input and, in 
addition, the simple and transparent form of this relation. Both features will 
be of great advantage in trying to pin down any possible inconsistencies of the 
Standard Model. A window for New Physics will exist above the largest possible 
Standard Model value implied by ( P7| ) and the model-independent upper bound 
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B{Kl Tc^iyu) < 4:A ■ B{K+ tt+z/z/) < 6.1 ■ 10"^ (50) 

Recent discussions on specific scenarios that might lead to an enhancement of 
K — > TTvv branching ratios over their Standard Model expectations can be found 



in [|T^] - [^. As pointed out recently in p3|, the CP- violating ratio e' /e plays 
in spite of large theoretical uncertainties a significant role in bounding possible 
enhancements. 

Next, we would like to point out an independent upper limit on B{Kl — ^ 
7r°i/z/), which follows directly from the upper limit on \ Vub\, using rj < \Vub/Vcb\/ \- 
It reads (kl/A^^ = 6.78 ■ 10"^) 

B{Kl ^ vrVz7) < ^ \V^bVcbX{xt)f < 6.5 ■ 10"" (51) 



Although still affected by the theoretical uncertainties in extracting \Vub\, this 
very clear and direct bound will become increasingly useful as our understanding 
of \Vub\ improves. 

Finally, we remark that once Kl —>■ n^uu and — > n'^uu have been mea- 
sured, the direct extraction of sin 2/3 and ImA^ from these modes becomes possible 
and offers additional opportunities for accurate tests of flavour physics M . 
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7 Conclusions 



The Standard Model predicts a characteristic pattern of rare, flavour-changing 
neutral current processes. Different transitions are interrelated in particular 
through the underlying common CKM structure of quark mixing. To make de- 
cisive tests of the flavour sector it is mandatory to have flrm control over all 
theoretical uncertainties. 

In the present paper we have discussed specific examples of relations that rely 
only on theoretically clean observables. For instance, knowledge of AM^/ AM^ 
and sin 2/5 determines B{Kl — ^ n^vv). An upper limit on AM^/ AM^ alone yields 
a stringent upper limit on B{K^ — > tt+i/z?). These constraints are both simpler 
and theoretically cleaner than the results of a standard analysis, where CKM 
parameters are determined using Ek, \Vub/Vcb\ and AM^. Such an analysis is at 
present still very useful to obtain benchmark predictions for rare decays. We have 
updated previous analyses for K nuu |jT6|, obtaining 

^ TT+z/p) = (0.82 ± 0.32) ■ 10"^° 
B{Kl n'^iyu) = (3.1 ± 1.3) ■ 10"" 

However, in the future one might want to avoid the hadronic uncertainties neces- 
sarily affecting these estimates, using instead the relations mentioned above. We 
note especially that a very clean, and conservative, upper bound 

BiK+^-n+vv) < 1.67-10"^° 

can be deduced with present knowledge from AM^/AMs alone. Higher values 
would clearly be in conflict with the Standard Model. It will be interesting to 
follow the future development of the experimental result on B{K^ — > tt'^uu) in 

The clean constraints on B{K^ —>■ -k^vv) and B{Kl —>■ ti^vv) we have dis- 
cussed should be useful to guide the interpretation of future experimental data 
and might well be essential to establish, with confldence, the existence of New 
Physics. 
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Appendix 



Top Contribution 

In the following we discuss the origin of the additional term that needs to be 
included in our previous results in 0. It arises from the box-diagram contribu- 
tion to the effective hamiltonian and is related to the existence of an evanescent 
operator (non-vanishing only in Z) 7^ 4 dimensions) within the framework of di- 
mensional regularization, as pointed out in . In the following argument we will 
use a gluon mass to regulate infrared (IR) divergences. A different IR regulator 
(external quark masses) was used in 0. Both approaches lead to the same result. 

To be specific we first concentrate on the top contribution for the rare decays 
with neutrinos in the final state. If we define 

Q = sYl'^Yil - 75)rf^>7^7a7.(l - 15)1^ (A.l) 

Q = S7^(l - 75)rfz77^(l - 75)z/ (A.2) 

we may write 

Q = {lQ + ae)Q + E (A.3) 

which defines the evanescent operator E. Here e = {A — D)/2 and a is an arbi- 
trary constant number. To lowest order in the one-loop box graph is simply 
proportional to Q. It therefore contributes to the hamiltonian (in D dimensions) 

Hto. ~ ^BoQ = BQ + BeE (A.4) 

where 

B = B^ + 0{as, e) Be = 4^0 + (A.5) 

lb 

Writing B^ = {B,Be), = {Q,E), and including ultraviolet (UV) counter- 
terms, ( |A.4| ) becomes 

Hfeox ~ B^Q + B^iZ2Z-' - 1)Q (A.b) 

Here Z2 is the quark-field renormalization constant and Z the renormalization 
constant matrix of Q defined by 

(Q)W = Z,'Z{Q) (A.7) 

relating unrenormalized (*^°^) and renormalized amputated Green functions with 
operator Q insertion. 

To 0{as) the matrix Z can be calculated by considering one-gluon exchange 
across the quark current in Q and E. To extract the UV renormalization, we may 
for instance treat the external quarks as massless and on-shell and use a gluon 
mass A as infrared regulator. In this case 
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and the matrix Z is found to be 




r = 32 (A.9) 

Z implies a finite renormalization of {Q,E). Since the quark component in Q 
is a conserved current, Q by itself is not renormalized in ( [A.9| ). On the other 



hand, E mixes back into Q under finite renormalization ~ r. Note that r is 
a pure short- distance quantity. It arises when the ultraviolet pole of the loop 
integration ~ 1/e multiplies a term ~ eQ from the D-dimensional Dirac algebra. 
Note also that the arbitrary number a in ( |A.3| ) does not appear in relation ( [A.9|) 



and is unimportant in the following. The element s in ( [A.9| ) is likewise irrelevant 



in the present context. The hamiltonian in (|A.6|) may then be rewritten as 

Hto. ~ (5 - ^tBe^ Z2Q + BeZ^E + . . . (A.IO) 

where irrelevant terms have been omitted. 

To extract the Wilson coefficient 5, the matrix element of the effective theory 
( |A.10 ) has to be matched to the corresponding (renormalized) amplitude in the 



full theory, Z2F. If we denote the tree-level matrix element of operator O as {0)t 
and use a gluon mass A as infrared regulator, we have 

Z2{QY'^ = {Q)t Z,{E)^'^ = ^r{Q)T + . . . (A.ll) 

and 

Z2F^''\Q)t = B{Q)t (A.12) 



In this case the contribution of r has canceled on the r.h.s. of ( |A.10| ) and the 



coefficient B is given directly from the calculation in the full theory, B = Z2F^^^ 
We see that the presence of the evanescent operator is then irrelevant and can be 
ignored from the beginning. This is not unexpected, because the l.h.s. of ([A.12|) 



is finite after renormalization and the limit D — 4 can be taken without the need 
to specify the subtraction of evanescent operators (in contrast to the case of the 
ordinary non-leptonic hamiltonian |^). This situation was implicitly assumed 



in the calculation of the 0{a^ box contribution to the Wilson coefficient in [0. 
However, in IR divergences were regulated dimensionally, rather than by a 
gluon mass. This is a legitimate procedure, but it involves a subtlety that was 
not taken into account in the original work 0|. 

In the case of dimensional regularization of both UV and IR divergences, and 
with external quarks massless and on-shell, Z2 = 1 and the 0{a^ corrections to 
the matrix elements of the operators in ( |A.10| ) vanish identically. The matching 
relation between full and effective theory ( [A.10| ) becomes, instead of ( |A.12| ) 



F(^)(g)T= (5-^r5^)(g)T (A.13) 
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The relevant coefficient B is therefore not given by the full theory amplitude alone 
(unlike the case of (^D), but by F^^^ + as/{^T^)rBE. With the QCD corrected 
box function written as p| 



B{x, 1/2) = 5o(x) + 1/2) (A.14) 

in the case of external neutrinos (T3 = 1/2), this implies that an additional 
term Ai?i(l/2) = tBe = 2i?o has to be added to the result in 0. Recalling 
X = C — 45(1/2), this corresponds to AXi = — 8i?o, which is already included 
in (H). It can be checked that the correction term as/ {^T^)rBE is identical to 
the (finite) difference between the full theory amplitudes with gluon mass and 
dimensional regularization, respectively, of IR sing ularities, i.e. ZgF^^) - F^^). 
Thus, identical results are obtained for B independent of the IR regulator, as it 
must be the case. 



Charm Contribution 

A similar correction arises also in the charm sector. To determine the correspond- 
ing change in the box function we write the result of the renormalization group 
analysis as 



(1/2) _ x{m) II 



m? 1 - r 256 6 625 I 3750 256 

Here we have explicitly displayed the dependence on 7^2^ . The element 712'* of the 
two-loop anomalous dimension matrix for the box contribution determines the 
NLO correction in B^l (apart from the anomalous dimensions of and of the 
charm quark mass, which are well known; more details on the renormalization 
group analysis can be found in 0). 7^2^ can be fixed by expanding B^l to ffist 
order in and comparing with the top contribution B{x) expanded to ffist order 
in X. This yields 

7ff = f (A.16) 



and one obtains (|T0|) from ( |A.15| ) 



We next show that the same result follows also directly from an analysis 
of operator mixing. This discussion extends the NLO calculations for the box 
contribution in the charm sector of ir'^uu, presented in 0, to include the 

effect of the evanescent operator E ( |A.3| ). 

The relevant basis of operators is given by O = (O, Qc, E^)'^ with 

= -ij d^xT{{sc)v-A{yl)v-A) {x) {{h)v^A{cd)v-A) (0) -{c^u} (A.17) 
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2 2 
771 Til 

= A^Q E, = ArE (A.18) 

where m = rfi(.{^) and g is the QCD couphng. 

The ingredients of the renormahzation group analysis are the initial values 
of the Wilson coefficients and the anomalous dimensions governing the evolution 
from scale /i = Mw to fi = rric- Details on the procedure may be found in |^ 0. 
Here we concentrate on the role played by the evanescent operator. 

The initial conditions are not affected by the presence of Ef.. The coefficient 
of Ec itself is 0{as), which gives an irrelevant contribution at NLO (where the 
highest order terms kept are 0{as) times a physical operator). However, E^. con- 
tributes to the physical (O, Qc) sector of the anomalous dimension matrix through 
operator mixing, affecting the NLO quantity 712'' • We work with dimensional reg- 
ularization throughout, for both UV and IR divergences. The renormahzation in 
the {Qci Ec) sector can be written as 

^2(gc)^°) = {Qc) = {Qc)t (A.19) 

Z2{E,f^ = Zs2{Qc) + ZssiE,) = {E,)t (A.20) 

Unrenormalized, renormalized and tree-level matrix elements are denoted by {)^'^\ 
and ()t, respectively. Here and in the following the renormahzation of m and g 
is, however, understood to be already carried out in the matrix elements labeled 
by 



From ( |A.9| ) above we have 

Z32 = ^r Z33 = 1 + (A.21) 

Next 

Z2(0)(°) = (O) + Z^2{Qc)t + Z^3{E,)T (A.22) 



where 



(O) + {Zu + ^13^32) (Qc) + ^13^33 (^c) 



Zu = f^ + 0{al) Z,s = l^ + 0{al) (A.23) 



The renormahzation of O = (O, Qc, E^.)'^ may be summarized as 

Z2(d)(°) = (1 + 0(0) (A.24) 

/ Zi2 + Z13Z32 Z13Z33 



C 



I = 0(a,) (A.25) 

\ Z32 ^33 ~ 1 



The anomalous dimension of O is given by 



7 = (1 + C)-'r(i + c) + (1 + C)-':^ (A.26) 
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where F = diag(0, 722, 722), with 722 the anomalous dimension of m?/g'^. Using 



d 



dliafi 



Z2(d)W = -rz2(d)(°), 



d 



(iln/x 



(d) = -7(d) 



(A.27) 



( [A.26| ) follows from ( [A.24| ). To see the impact of Ec (entries with index '3') on 
7i2, we note that the quantities 722, -^'is and Z32 all start at 0{as). We then 
find 



. , dCi2 . dC32 
7l2 = -C12722 + -T, Cl3 

ct m/i 

= -Z12722 + — h 



dlinjj dlnjj 
The contribution from the evanescent operator is 

dZvi ^ 



dliifx 
^^^^ Z32 + 0{al) 



As7i2 



c/ln/i 



^32 



-2r 







(sr-« 






KAtt) 



where we have used ( |A.21 ), ( A.23| ) and das/dlnfi = —2eas + 0{al 
two terms in ( [A. 291 ) were computed in With 



7l2 



Ois (0) , 



An) 



we then obtain, using the result of 



320 . (1) 



128 



(A.28) 
(A.29) 

(A.30) 

The first 

(A.31) 
(A.32) 



in agreement with (|A.16| ). 



Summary 

In the two preceding subsections we have discussed the case with final state 
neutrinos (T3 = 1/2, X-functions). Completely analogous arguments apply to 
the case of external charged leptons (T3 = —1/2, y-functions) and lead to the 
results summarized in the main text. We finally collect, for the sake of clarity, all 
changes due to the additional contribution from the evanescent operator. They 
affect only the box contribution. We list the terms that have to be added to 
our previous results in [g, ^, |^ for the corresponding Wilson coefficient functions. 
The modifications in the top sector read explicitly 



AX{x) = -8Bo{x[ 



as 
An 



AY{x) = +8Bo{x 



An 



(A.33) 



where Bq is the leading order box function 



Bn(x) 



X xlnx 
+ v. 



X 



x] 



(A.34) 
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In the charm sector one has 

1 24 1 24 

AXnl = ~-xim)Kc'' {KK2 - 1) AYnl = +-x{m)Kc'' {KK2 - 1) 

(A.35) 

in the notation defined in sections 2 and 3 of the main text. For both Xj^^ and 
Yatl the relevant entry of the two-loop anomalous dimension matrix 7^2'' changes 
by A^'S = -64. 
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